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ANALYSIS OF HELICOPTER MANEUVER-LOADS 


AND ROTOR-LOADS FLIGHT-TEST DATA 

By Edward A. Beno 
Sikorsky Aircraft Division 
United Aircraft Corporation 


SUMMARY 

A study was conducted in which available airload and blade response data 
for the NH-3A and CH-53A rotors were analyzed in an attempt to provide greater 
insight into the sources of rotor vibratory loads in both level and maneuvering 
flight. Primary emphasis in the study was placed on examining and understand- 
ing causes of high-frequency rotor control loads. Secondary objectives were 
(l) to examine the effect of number of rotor blades on hub vibratory shear 
forces and (2) to assess which of the many terms appearing in the hub vibra- 
tory shear force expressions were of most significance. 

Although stall flutter is probably a major contributor to high-frequency 
control load oscillations when operating well above the knee of the control 
load curve, it was found for the less stalled conditions analyzed in this study 
that the principal high-frequency control load cause appears to be a forced dy- 
namic response of the blade to the time varying aerodynamic moment. Other fac- 
tors previously considered to be potentially important, such as moments induced 
by blade/tip-vortex interactions or blade bending, were generally small for the 
conditions studied. 

Analysis of the hub vibratory shear forces indicated that the nondimen- 
sionalized force harmonics for the six-bladed rotor were generally higher than 
those for the five-bladed rotor. Care must be exercised in generalizing this 
result, however, inasmuch as detailed blade dynamic characteristics (e.g. 
natural frequencies) undoubtedly influence the comparison. Finally, a review 
of the equations used to compute hub vibratory forces from measurements of 
blade root dynamic response indicates that many of the dynamic correction terms 
due to motion of the hub, blade cuff, and instrumented blade segment can be 
neglected. 


INTRODUCTION 


One of the current problems facing helicopter rotor designers is the 
flight envelope limitations imposed by high-speed and maneuvering flight con- 
ditions where high vibratory control loads are encountered. Expansion of the 
envelope may be possible if the source or sources of these vibratory loads can 
be established and accurate methods validated for their prediction. 

Previous investigators have studied this problem including analysis of 
the effects of stall flutter, rotor inflow, and unsteady aerodynamics. The 


present study sought an explanation of the cause of the vibratory control loads 
measured on the NH-3A and the CH-53A, and the ability of available analytical 
methods to predict the measured blade response. References 1 and 2 have con- 
cluded that retreating blade oscillations can be the result of a complicated 
series of dynamic and aerodynamic events. Reference 3 postulated the source of 
high oscillatory control loads to be a wake-induced forced response of the tor- 
sional degree-of- freedom. Vortex crossings in the wake were shown to have a 
direct correspondence with measured retreating blade c m oscillations. The study 
reported herein investigated torsional response utilizing in-flight data pre- 
viously measured on the NH-3A and the CH-53A. Reference describes the measure- 
ment program on the NH-3A. The CH-53A data were acquired during a previous flight 
test program funded by the Naval Air Systems Command under Contract NOw 63-0150-f, 
CO-kO (63-0150) and CO-^ORl (63-0150). 

Included in the present study was an analysis of the measured maneuvering 
and high-speed blade azimuthal responses. Analysis was limited to a relatively 
few rotor revolutions of flight data. An understanding of the source of re- 
treating blade torsional oscillations was sought and measured blade response 
was correlated with analytical calculations. 

A secondary investigation covers an examination of the vibratory rotor 
loads on both the NH-3A and the CH-53A. Methods were sought from which the 
rotating blade hinge force harmonics and the non-rotating vibratory hub forces 
could be approximated with good accuracy from measurements at the root of one 
blade. Comparisons of the hinge force harmonics were made between two rotor 
systems having a different number of blades. 

Robert H. Blackwell, Jr., Sikorsky Aircraft Division, United Aircraft 
Corporation, assisted in performing the correlation of test data with analysis. 


SYMBOLS 

Values are given in both the International System of Units (Si) and U. S. 
Customary Units. The measurements and calculations were made in U. S. Customary 
Units . 

b number of main rotor blades 

c blade chord, m (in.) 

c m section pitching moment coefficient, m 

% pc z [ftR(r/R + ysini/j )] 2 

c n section loading coefficient, l 

2§pc[fffi(r/R + ysini /; )] 1 

C t/o rotor thrust coefficient/solidity ratio 

C.I. constant inflow analysis 
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blade offset, m (in.) 
longitudinal rotor hub force, N (lb) 
radial blade hinge force, N (lb) 

blade radial force at inboard end of instrumented segment, N (lb) 

lateral rotor hub force, N (lb) 

edgewise blade hinge force, N (lb) 

vertical rotor hub force, N (lb) 

vertical blade hinge force, N (lb) 

frequency relative to cutoff 

Gaussian coefficients 

mass moment of inertia of instrumented segment about flapping axis 
N-m-sec2 ( lb-in. -sec 2 ) 

mass moment of inertia of instrumented segment about lag axis, 
N-m-sec 2 (lb-in. -sec 2 ) 

blade torsional mass moment of inertia, N-m-sec 2 ( lb-in .i-sec 2 ) 

section aerodynamic pitching moment factors, nr (in. ) 

blade root torsional spring constant, N-m/rad (in. -lb/rad) 

section normal loading, N/m (lb/in.) 

root torsional structural moment, N-m (in. -lb) 

aerodynamic pitching moment, N-m (in. -lb) 

model root torsional structural moment, N-m (in. -lb) 

blade flatwise bending moment, N-m (in. -lb) 

blade edgewise bending moment, N-m (in. -lb) 

advancing blade tip Mach number (r/R = 1.0, \p = 90 deg) 

section aerodynamic pitching moment, N-m/m (in. -lb/in.) 

blade mass inboard of instrumented segment, kg (lb-sec 2 /in. ) 

blade instrumented segment mass, kg (lb-sec 2 /ln. ) 
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n/rev 

P 

P ~,h 

R 

r 

T 

V 

V.I. 

XI, X2, . . 


frequency of n-th harmonic, i.e., n times the main rotor frequency 
absolute pressure, N/m^ (psia) 

p 

free stream static pressure, N/m (psia) 

rotor radius, m (in.) 

blade radial location, m (in.) 

modal contribution of aerodynamic pitching moment, N-m (in. -lb) 
forward airspeed, m/sec (kn) 
variable inflow analysis 

symbolic blade radial hinge force terms (see table II ) 

radial coordinate in blade axis system of center of gravity of 
instrumented segment, m (in.) 

radial distance, in blade axis system, from inboard end of instru- 
mented segment to center of gravity of instrumented segment, m (in.) 


x 

c 


radial coordinate, in blade axis system, of center of gravity of 
mass inboard of instrumented segment, m (in.) 


Xj^ rotor head displacement in shaft rotating system radial direction, 

m (in.) 

Yl, Y2,.. symbolic blade edgewise hinge force terms (see table II) 


y 


blade chordwise location measured from leading edge, m (in.) 


edgewise coordinate, in blade axis system, of center of gravity of 
mass inboard of instrumented segment, m (in.) 


'R 


rotor head in-plane displacement in shaft rotating system perpendi- 
cular to hub radial direction, m (in.) 


Z1,Z2, . . 


symbolic blade vertical hinge force terms (see table II ) 


z 

c 


vertical coordinate, in blade axis system, of center of gravity of 
mass inboard of instrumented segment, m (in.) 


J R 


rotor head vertical displacement in shaft rotating system, m (in.) 


e 


blade root flapping angle, deg 


Y 


blade root lag angle , deg 
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Ap 

Ax 

S 

0 

0 T 

V 

P 

Y n 

ip 

n 

Subscripts : 


differential pressure, N/m^ (psi) 

radial length of blade instrumented segment, m (in.) 
damping ratio 

blade root pitch angle, deg 

blade root torsional elastic twist, deg 

advance ratio, _V 

mass density of air, kg/m ^ (lb-sec^/in. ) 

phase angle associated with n-th harmonic, deg 

blade azimuth angle, deg 

rotor angular velocity, rad/sec 

blade torsional natural frequency, rad/sec 


( ) c cosine component of harmonic 

( ) MT refers to summation of major terms 

( ) n-th harmonic order 

n 

( ) g sine component of harmonic 

( ) 1 at the inboard end of the blade instrumented segment 

( )^ at the outboard end of the blade instrumented segment 

) refers to 5-bladed NH-3A rotor system 

g( ) refers to 6-bladed CH-53A rotor system 

Miscellaneous : 

( ) first derivative with respect to time 

( ) second derivative with respect to time 


DATA ACQUISITION AND REDUCTION 


The data analyzed in the study reported herein were measured 


vious NH-3A and CH-53A flight test programs. The flight tests were primarily 
performed to acquire data describing the aerodynamic and dynamic characteris- 
tics of each aircraft. 

The NH-3A is a compound helicopter capable of not only high forward speeds 
but also of varying rotor lift and drag. Two auxiliary turbojet engines and a 
wing including a full-span flap provide the aircraft with these capabilities . 

The 18.90: m (62 ft) diameter main rotor is a 5-bladed fully articulated system. 
The -1* deg twisted blades have a modified NACA 0012 airfoil section. Normal 
gross weight is 8560 kg (189OO lb). A photograph of the NH-3A is included as 
figure 1. 

The CH-53A, a high speed assault transport helicopter, is capable of 
flight at higher gross weights than is the NH-3A. The aircraft was flown at 
both 15860 kg (35OOO lb) and 19040 kg (42000 lb)’ gross weights. The 6-bladed 
fully articulated CH-53A main rotor has a diameter of 21.95 m (72 ft) and the 
-6 deg twisted blades have a modified NACA 0011 airfoil section. The CH-53A is 
shown in figure 2. 


Instrumentation 

An airborne data acquisition system was used to record the aerodynamic 
and dynamic characteristics of each aircraft. Approximately one hundred chan- 
nels of data were recorded on a magnetic tape recording system installed in the 
cabin of the aircraft. This harrow band system provided a continuous frequency 
modulation multiplexed recording of each measurement. Signal conditioning of 
the multiplexed output signal was performed before the data were actually re- 
corded on the tape transport. Quasi-static data describing the recorded flight 
condition were provided on 35 mm film that recorded photopanel measurements . 

The photopanel measurements included airspeed, altitude, rate of climb, rotor 
angular velocity, ambient atmospheric conditions, and aircraft attitude. A 
complete description of the instrumentation and data acquisition methods used 
is included in reference 4 . 

A dynamic data sampling signal, used for time reference, was provided in 
the aircraft by a rotor azimuth pulse generator which divided a main rotor rev- 
olution into 5 deg increments. A slotted disc and optical system mounted on 
the bottom of the main rotor shaft provided this 72/rev signal plus a l/rev 
pulse used as the 0 deg azimuth reference. This corresponds to the instrumented 
main rotor blade aligned over the aircraft tail. 

A block diagram of the airborne recording system used on the CH-53A is 
shown in figure 3. Basically the same system was utilized to record the NH-3A 
data. The measured blade response data used in the study reported herein in- 
clude blade pressures; flatwise, edgewise, and torsional bending moments; and 
pushrod loads. The distribution of the miniature absolute and differential 
pressure transducers on the instrumented blade is given in figure 4 . Absolute 
pressure transducers were used on the upper and lower surfaces of the blade 
spar while differential gages were installed through the blade pockets and tip 
cap. Bending moment measurements were distributed along the blade span to re- 
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cord the response of the blade. 


Data Processing 

The primary objective of this program was to perform a detailed analysis 
of a limited number of rotor revolutions of important flight data. The emphasis 
was on high-speed and maneuvering flight conditions for both the NH-3A and the 
CH- 53 A. The maneuvering conditions considered in this program were comprised 
of banked right-turn regimes. Left turn data were not available. For the 
WH-3A, the measured data were acquired only during the actual steady-state 
period of each banked turn. The available CH-53A right-turn conditions, how- 
ever, included data measured in the trimmed level flight regime, through a 
transition period, and finally during the actual steady-state portion of the 
coordinated turn. Therefore, the NH-3A data cycle studied for the turns con- 
sisted only of a rotor revolution during the actual banked turn. For the CH-53A 
maneuver, three separate cycles were investigated; trimmed level flight, tran- 
sition, and steady-state coordinated banked turn. Table I describes each of 
the data cycles investigated in the high -frequency blade response part of this 
study. Case numbers refer to those used in reference k and earlier CH-53A 
studies. ' 

Processing of the measured data was initiated by the formation of oscil- 
lograph traces from the flight test magnetic (analog) tapes. These records 
were arranged to include the 1/rev and 72/rev signals required for subsequent 
data reduction. In order to retain possible high-frequency response of the 
test data, a filtering cut-off frequency of 100 Hz was utilized during the gen- 
eration of the oscillograph rolls. This had the effect of filtering out un- 
wanted signal noise above 100 Hz while retaining real signal response at fre- 
quencies up to approximately 30/rev, depending on the rotor angular velocity. 
The signal amplitude response using the low pass filters having a cut-off fre- 
quency of 100 Hz is shown in figure 5. 

Selection of the particular cycles of data to be used in the high-fre- 
quency blade response study were chosen by manual scanning of the oscillograph 
traces. For maneuvers and high-speed conditions the cycles chosen within each 
regime, described in table I, were selected primarily on the basis of their 
high-frequency torsional response content. For these cycles the oscillograph 
roll data for each required measurement were reduced to punch card form. 

The 72/rev signal included on the traces was used as a time step interval 
at which the data were punched into cards. Each parameter was punched as a 73 
point time history corresponding to a 0 to 360 deg rotor azimuth period ( 5 deg 
azimuth interval). Along with the actual dynamic data for each measurement, 
calibration factors from the traces were also punched on the cards for later 
use in conversion of the machine units to engineering units. The dynamic data 
and calibration factors on the punch cards were input to a computer program for 
final processing. 

A secondary study included in this report was an analysis of main rotor 
loads. Basic objectives were to determine the major dynamic terms that contri- 
bute to harmonics of blade hinge forces and b/rev hub loads and the effect of 
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the number of main rotor blades on harmonics of the hinge forces. This was ac- 
complished by an analysis of the existing data measured during the flight pro- 
grams on the NH-3A and CH-53A. Cases were chosen from the available data to 
cover a wide range of operating conditions. 


Data Analysis Methods 

Digital computer programs were utilized for final processing of the 
measured flight-test data. The blade pressure, bending moment, and pushrod 
load data punched on cards from the oscillographs were converted to engineer- 
ing units and listed by one program. A second digital program was used during 
the original flight-test programs to calculate the blade hinge forces and non- 
rotating b/rev hub loads. From the measured data, the separate dynamic terms 
contributing to these loads could be calculated by an option in this program. 

With laboratory and aircraft calibration factors, along with the oscillo- 
graph conversion factors punched on the cards, the machine unit blade response 
data on punched cards were first converted to engineering unit data by a com- 
puter program. This procedure gave the azimuthal history of the measured pa- 
rameters including response at frequencies up to 100 Hz (filtering cut-off fre- 
quency). For each aircraft, 33 measurements were listed in this form. These 
included 2k blade absolute and differential pressures at 75 , 85, and 95 percent 
radial stations; 3 flatwise bending moments; 3 edgewise bending moments; 2 tor- 
sional moments; and the pushrod load. At stations where two absolute pressure 
transducers were mounted, on the upper and lower blade surfaces, the differen- 
tial pressure time histories were also calculated. At each of the three blade 
radial stations the measured chordwise distribution of pressure was then uti- 
lized to calculate the blade section loading (airload) and aerodynamic pitching 
moment time histories. The five chordwise pressure stations were at 4.2, 15.8, 
30.0, 60.0, and 91.0 percent of chord. A Gaussian integration scheme was used 
to calculate the time history of blade airload from the relation 


Z = c I (Ap). , G. 
* i=l 1 »* 1 


(1) 


where G. are the Gaussian coefficients for the five chordwise stations from G, 

l 1 

corresponding to k.2 percent chord to G<- at 91*0 percent chord. The coefficients 

are G.^ = 1/10, G^ = 1/10, G^ = 2/9, G^ = 16/45, and G^ = 2/9. The aerodynamic 

pitching moment about the blade feathering axis was also calculated from the 
pressures utilizing a trapezoidal distribution between the measured pressures. 
Leading- and trailing-edge differential pressures were assumed to be zero. The 
relationship for calculation of the pitching moment is 


m 


* 


5 

E 

i=l 


< 4p > i,* K i 


( 2 ) 


where Kj_ are blade geometry factors corresponding to the five chordwise stations 
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from at 4.2 percent chord to at 91-0 percent chord. For each aircraft the 
factors are 


Factor 

NH-3A 

CH-53A 

2 

m 

. 2 
m. 

2 

m 

. 2 
in. 

K i 

3.11 X10 -3 

1.82 

6.32 X10' 3 

9-79 

K 2 

2.31 

3.58 

1.69 

7.27 

K 3 

-1.88 

-7.56 

-9.90 

-15 . 34 

K 4 

-23.16 

-35.89 

-47-00 

-72.85 

K 5 

-25.21 

-39-08 

-51.17 

-79.32 


From the section loading and aerodynamic pitching moments calculated in equa- 
tions (l) and (2), the azimuthal histories of section loading and pitching mo- 
ment coefficients were calculated. 

The digital computer program used to calculate the dynamic terms contri- 
buting to the blade hinge forces and b/rev hub loads is described in detail in 
reference 4. The Root Response Method is programmed to calculate azimuthal 
histories of the blade hinge forces from measurements made at the root of the 
blade. These include flatwise and edgewise bending moments, blade tension, 
blade root angles, and accelerations at the blade hinge. The resulting equa- 
tions, given in reference 4 include a large number of dynamic terms in each of 

the three orthogonal hinge forces: radial, F*; edgewise, F*; and 

x y 

vertical, F*. These terms are harmonically analyzed. Assuming each main rotor 
blade encounters the same loads azimuthally as the instrumented blade, the non- 
rotating system b/rev loads are then calculated. This computer program provided 
the calculation of the separate terms required to determine the major terms con- 
tributing to the hinge forces and non-rotating loads. This program also pro- 
vided comparisons of the effects of the number of blades on these loads for both 
the NH-3A and CH-53A measured data. 


ANALYSIS OF HIGH-FREQUENCY BLADE RESPONSE 


The output of the digital computer program used to process the test data 
read from oscillograph traces provided the basis for a timewise study of high- 
frequency blade torsional response. It also made possible comparisons of dif- 
ferent conditions and different aircraft. Azimuthal time history plots referred 
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to in subsequent sections of the report are of two forms. Plots of level-flight 
conditions given in table I indicate an azimuthal range from 0 to U00 deg. For 
these plots the first 40 deg of the rotor revolution have been repeated to en- 
able better visualization of any high-frequency blade response oscillations oc- 
curring near the tail region of the aircraft. Maneuvering conditions are plotted 
only from 0 to 360 deg since for these regimes, the 0 and 360 deg azimuth 
values are generally different. 


Factors Affecting Blade Torsional Response 

The occurrence and magnitude of high-frequency blade torsional oscillations 
was a primary objective of this study. To determine the effect of the operating 
condition on this response, various conditions were compared. Of specific in- 
terest was the buildup of loads either at high airspeeds or during maneuvering 
conditions. Analysis of the measured blade response included studies of blade 
pressures, section loading coefficients, section pitching moment coefficients, 
blade/tip-vortex wake interactions, blade bending, and aerodynamic pitching 
moments. The effects of these parameters on blade torsional response are dis- 
cussed in the following sections. 

Effects of airspeed. - The effects of airspeed on the response of the NH- 
3A blade are shown in figure 6. On the advancing side of the rotor disc the 
section loading coefficient, c n ; section pitching moment coefficient, c m ; and 
the root torsional moment, M x , are basically unaffected by a change in airspeed. 
On the retreating side, however, moderately higher c n and c m amplitudes general- 
ly exist at the higher airspeed but the frequency content is the same. The 
change in is considerably different than that in c n and c m due to increased 
airspeed. While the steady value of the root torsional moment remains approxi- 
mately the same , there is a large increase in the vibratory content . The fre- 
quency of the oscillations is between 6/rev and 7/rev which corresponds to the 
blade torsional natural frequency of 6.6/rev for the operating conditions shown 
in figure 6. The measured pushrod loads for these conditions similarly indicate 
a buildup of the retreating blade oscillations at the same frequencies. 

Effects of maneuvers . - To determine the effects of maneuvering flight on 
blade response, CH-53A cases UlLF (level flight), i+1 TR (transition into a 
banked turn), and Ul TN (steady-state banked turn at a 60 deg angle of bank) 
were investigated. These effects are shown in figures 7, 8, and 9. An \increase 
in the steady value of c n can be noted in the typical data, shown at r/R^= 0.85 
in figure 7, as the aircraft progresses from level flight, through a transition 
period, and finally into the 2g coordinated turn. This corresponds to increased 
rotor thrust required to maintain attitude. Also shown is an increase in high- 
frequency content in the first and fourth quadrants of rotor azimuth for the 
turn. This is due to an increased effect of the instrumented blade passing 
closer to the trailing vortices shed from preceding blades. Larger magnitude 
oscillations occur in the fourth quadrant compared to those in the first. 

Typical of the spanwise distributions for the same CH-53A maneuvering con- 
ditions, are the section pitching moment coefficients at r/R = 0.85 given in 
figure 8. Compared to level flight, the steady values of c m are not appreciably 
affected by the maneuvering conditions, but a definite increase in the higher 
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harmonic content is evident. The advancing side of the rotor disc experiences 
several oscillations of very high frequency, above 30/rev, which also are prob- 
ably due to blade vortex crossings in the wake. In the vicinity of 270 deg az- 
imuth there is a growth of one small amplitude oscillation during level flight, 
into a short duration high amplitude spike during transition, and finally into 
two separate high magnitude fluctuations in the turn cycle. Along with the gen- 
eration of the two spikes near = 2J0 deg, a third large amplitude fluctuation 
appeared in the region over the aircraft tail during the turn. Although instru- 
mentation attitude limitations prevented knowledge of the exact rotor orienta- 
tion with respect to the free stream during turns, it is believed that the large 
amplitude c m variations are due to vortex crossing effects. Examination of the 
measured blade pressures at r/R = 0.85 indicates that large excursions in the 
differential pressure at the 4.2 percent chord station were the cause of the 
corresponding c m fluctuations. These are discussed in a subsequent section of 
this report. 

The c m data shown in figures 6 and 8 have indicated a steady leading edge 
up condition existing for both the NH-3A and CH-53A conditions. Steady c m values 
calculated from measured pressures vary between 0.02 and 0.04. Existing stan- 
dard NACA 0011 and 0012 steady and unsteady airfoil data, references 5 and 6, 
indicate that these positive steady values are unattainable over the angle of 
attack and Mach-number environment experienced in the flight conditions. The 
basic NH-3A blade has a modified 0012 airfoil and a modified 0011 airfoil is 
used on standard CH-53A blades. However, both of the test instrumented blades 
had chordwise epoxy fairing bands covering the area around the surface mounted 
pressure transducers and a spanwise fairing over the instrumentation wiring to 
the root of the blades. Details of the application and extent of these fairings 
are given in reference 4. It is possible, therefore, that these fairings, which 
had the effect of making the airfoil section unsymmetrical, may be the source of 
resulting leading edge up pitching moments on the blade. 

The existence of the positive pitching moments is also indicated in the 
measured blade root torsional moments shown in figure 9- Increased harmonic 
content is also evident in the torsional moment as the steady-state turn cycle 
is attained. The high amplitude retreating blade torsional moment oscillations 
during the turn do not indicate a direct correlation with the c m fluctuations 
shown in figure 8 due to variations of the c m distribution along the blade. 

The root torsional moment measurement records the integrated effect of loadings 
along the entire blade length. 

Effects of blade /tip-vortex wake interactions . - For 4-blade d H-34 maneu- 
vering conditions, reference 3 indicated that a direct correlation existed be- 
tween high-frequency retreating blade oscillations and the passage of each blade 
by the tip vortices generated by preceding blades. In the present study, the 
available cases were investigated to determine whether this interaction existed 
for the 5-bladed NH-3A or the 6-bladed CH-53A. It should be noted that the num- 
ber of blades and the advance ratio can affect the frequency and strength of the 
blade loads induced by the blade/vortex interactions. Also the sensivity of a 
given blade design to such periodic excitation will, of course, depend on the 
blade natural frequencies. 
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The non-distorted helical pattern of the blade tip vortex trajectories is 
shown in figure 10 for NH-3A case 72. Shown is a view from above and perpendi- 
cular to the rotor tip path plane. Indicated in the figure are the plan view 
intersections of the 95 percent blade station with the wake pattern formed from 
the tip vortex of each blade. Three of these crossings occur in the fourth 
rotor quadrant and six in the first quadrant. These intersections represent 
possible blade/wake interactions, but the proximity of the blade and the tip 
vortices, along with the circulation strength, depend on the rotor orientation 
with respect to the free stream and the wake inflow velocities. 

The azimuthal histories of blade differential pressures at 75, 85, and 95 
percent radial stations are shown in figures 11, 12, and 13 for NH-3A case 72. 
Included within the plots are the azimuthal locations of possible blade/wake 
interactions determined from the helical wake pattern shown in figure 10. In- 
dicated are the intersections of the tip vortices with the blade radial station 
corresponding to the plotted blade pressures. These locations also indicate 
the directional sense of the tip vortex tangential velocity as the instrumented 
blade approaches. Figures 11, 12, and 13 do not indicate, for NH-3A case 72, 
any strong correlation between the possible first quadrant vortex crossings and 
the pressure response. In general the pressure variations across the chord ap- 
pear to behave similarly except at the 30.0 percent, and to a lesser extent at 
the 15.8 percent, chord station on the advancing blade at 85 and 95 percent ra- 
dial stations. These fluctuations, in the general shape of a square wave, ap- 
pear to be due to a critical Mach-number effect which will be discussed in the 
succeeding report section. They are not believed to be the result of any blade/ 
wake intersections. On the retreating side of the disc there is evidence of 
probable vortex crossing effects, the strongest occurring at approximately 280 
deg azimuth where relatively rapid pressure fluctuations correspond to a possi- 
ble crossing at the same time. However, other lower frequency pressure oscilla- 
tions occur where vortex crossings do not appear to exist, such as those that 
initiate at the 85 and 95 percent blade stations over the nose of the aircraft. 
Vortex crossings are unlikely in this region since, as can be seen in figure 10, 
the blades are approaching undisturbed air. 

Based on the pressures shown in figures 11, 12, and 13, the NH-3A right 
turn at a 60 deg angle of bank does not appear to exhibit a strong correlation 
between the possible blade/vortex crossings and the measured pressure response. 
However, this maneuver condition did display a limited effect while the steady 
state level flight conditions that were analyzed did not exhibit any noticeable 
effect of the blade /vortex crossings. 

In figure lU, the c n and c m time histories, calculated from the measured 
pressures at r/R = 0.95, shown in figure 13, are plotted along with the tor- 
sional moment measured at the root of the blade. Included in the c m plot are 
the possible vortex crossings with the 95 percent blade radial station deter- 
mined from the plan view of the non-distorted tip vortex trajectories. No ap- 
preciable oscillations occur in c n or c m on the advancing blade other than one 
near ip = bo deg which is attributable to the critical Mach-number effect dis- 
cussed in the following section. On the retreating side, however, several os- 
cillations occur in both c Q and c m , two of which agree fairly well in phase 
with, vortex crossings at 2o7 and 317 deg azimuth. Comparison of the measured 
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root torsional moment in figure l4 with the azimuthal locations of possible 
blade/vortex crossings shows that the torsional moment oscillations initiate at 
180 deg azimuth, prior to any possible crossings. In fact the largest torsional 
moment oscillation occurs prior to the first helix crossing with the 95 percent 
span of the blade. The c m distributions at 75 and 85 percent blade span (not 
shown) similarly do not indicate blade/vortex crossings at those stations to be 
a direct cause of the root moment oscillations. 

Analysis of other NH-3A maneuver cases resulted in the same general con- 
clusions as described for case 72. A limited degree of correlation existed be- 
tween possible blade/vortex intersections and the pressure and c m response. 

None of the cases examined, however, indicated a direct correlation between the 
blade/wake interactions and the NH-3A retreating blade oscillations measured in 
both the root torsional moment and the pushrod load. 

The effect of the blade/wake intersections was also examined for CH-53A 
conditions. The measured pressure distributions at 75, 85 , and 95 percent blade 
radial stations for CH-53A case 4lTN, a 60 degree angle of bank right turn, are 
shown in figures 15, l6, and 17. The corresponding blade and non-distorted wake 
vortex crossings are also indicated. A very definite high-frequency response 
occurs in the advancing blade pressures, especially at 75 and 85 percent span 
stations. The r/R = 0.75 station indicates approximately 8 advancing blade os- 
cillations, appearing in all chordwise stations, occurring at the same time as 
the vortex crossings and at the same frequency, approximately 30/rev. In addi- 
tion several retreating blade oscillations occur roughly in the region of possi- 
ble crossings. 

Maneuver case 4lTN, however is atypical among the data for the CH-53A in 
so far as the effect of blade/wake interactions on blade pressures is concerned. 
Examination of the pressure time histories for other CH-53A cases, including the 
level flight and transition cycles prior to the case UlTN turn, indicates little 
or no fluctuations in the region of possible crossings. Therefore, although the 
exact rotor tip path plane orientation in the free stream is unknown for the 
turn, it does appear that the instrumented blade encountered strong vortex ef- 
fects as the blade passed close to, or through, the trailed vortices. The rotor 
orientation and coning angle of the other cases was apparently such that the 
trailed vortices did not pass close to the blade. Also, since the turn condition 
is generating higher rotor thrust than the level flight conditions examined, 
the relative strength of the vortices should be greater. 

The results of the pressure fluctuations on c n and c m at 75 percent blade 
span for CH-53A case 4lTN are shown in figure 18 along with the measured root 
torsional moment. The occurrence of the advancing blade pressure oscillations 
at all five chordwise stations results in similar excitation of c n and c m . The 
root torsional moment in this region also exhibits some high-frequency response 
at about 25/rev, slightly lower than the approximately 30/rev pressure oscilla- 
tion frequency at both 75 and 85 percent span. The frequency difference may be 
due to differently phased pressure oscillations over the remaining blade area. 

On the retreating side of the blade, fairly good agreement exists between the 
higher frequency pressure, c n , and c m fluctuations at 75 percent span. The c n 
distribution also has two relatively large amplitude oscillations near \p= 245 
and ip = 330 deg. Examination of the pressures in these regions, shown in figure 
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15, indicates a phase difference between the five pressures across the chord. 

This is similar to the findings of reference 3 during a maneuver condition and 
is probably indicative of local dynamic stall effects. The phasing differences 
result in the large Cjj oscillations. The measured root torsional moment for 
the CH-53A case i+lTN turn exhibits frequency response in the 20/rev to 30/rev 
range whereas the NH-3A case 72 turn shown in figure lH does not. This is at- 
tributable to the effect of the blade/wake interactions which are not evident 
in the NH-3A case. The large amplitude root torsional moment oscillations on 
the case 4lTN retreating blade are somewhat different from the c m oscillations 
at 75 percent span due to different azimuthal histories of c m at other radial 
stations. The c m at r/R = O.85 is shown in figure 8. 

In summary, blade/wake interaction effects appeared to occur strongly in 
only one case investigated, a CH-53A banked turn. For both the NH-3A and the 
CH-53A, maneuver conditions appeared to be affected more by these effects than 
did the level flight regimes. The level flight CH-53A and NH-3A cases exhibited 
little or no effect, apparently due to either an absence of close rotor blade 
proximity to the wake or low vortex strength. The CH-53A turn did, however, ex- 
perience relatively strong first and fourth rotor quadrant effects on measured 
blade pressures, c n , c m , and root torsional moment. However, these effects do 
not appear to be the direct cause of either NH-3A or CH-53A high amplitude re- 
treating blade control load oscillations. 

Effects of advancing blade critical Mach number. - Examination of the 
measured blade pressure data revealed a change in the azimuthal history that 
occurred in several cases in the region of the advancing blade. Reference 7 
describes a similar change in H-3^ pressure data which appears as an approximate 
square wave caused by two step changes in the differential pressure. This ef- 
fect is demonstrated by the NH-3A case 72 pressure shown in figures 11, 12, and 
13. For this case the occurrence of the step changes is most pronounced in the 
measured pressure at r/R = O.85 and y/c = 0.300 as shown in figure 12. Evidence 
of the effect is generally seen only in the 30.0 and 15.8 percent chord station 
pressures for all NH-3A and CH-53A cases. Figure 12 shows the sequence of pres- 
sure changes for the y/c = 0.300 station to be a step increase at ip = 50 deg 
followed by a step decrease to the original pressure level at \p = 115 deg. 

These square waves always appear roughly symmetric about an azimuth point be- 
tween 80 and 100 deg. This corresponds to the blade environment of highest 
local Mach number. 

Figure 19 shows this effect on the measured differential and absolute pres- 
sures for NH-3A case 72. The greatest effect at 30.0 percent chord is seen to 
occur in the 85 percent span station while a lesser magnitude effect is seen at 
r/R =0.95. No effect is apparent at 75 percent. The last three plots of fi- 
gure 19 present the corresponding measured absolute pressures on the upper and 
lower blade surfaces at the three spanwise stations. The absolute pressures 
shown have a zero reference when the aircraft was in a static ground condition. 
Due to flight at altitude, the free-stream static pressure changed and the abso- 
lute transducer readings included this delta static pressure change. The free- 
stream static pressure, p ro is indicated in the plots. As indicated, the 
square wave effect in the differential pressure is due to a corresponding square 
wave occurring in the upper surface absolute pressure. The largest step changes 


in absolute pressure occur at r/R = 0.85 corresponding to the large effect seen 
in the differential pressure. The Karman-Tsien relation, reference 5, was used 
to predict the absolute pressure required on the upper surface such that the lo- 
cal Mach number reaches unity. This critical pressure is dependent on the free- 
stream Mach number. In general this indicated that the absolute pressures 
reached in the area of the first step change correspond to pressure at or above 
the critical pressure. At 75 percent span, where no step changes in pressure 
occur, the Karman-Tsien relation predicts the critical pressures to be above 
the pre-sures on the blade. 

It is believed that the generation of the step changes in pressure are due 
to a shock -occurring over the upper blade surface when the critical Mach number 
is reached locally. The first step change in the negative absolute pressure 
sense could be due to the pressure buildup forward of the actual shock followed 
later in azimuth by a step change in the opposite pressure sense. The actual 
shock may be strongest aft of the 30.0 percent chord station. Due to the limit- 
ed number of pressure measurements across the chord, the next station was at 
60.0 percent, and the transducer at this station did not experience the step 
change effect. Analysis indicates that the effect does occur in the region of 
high free-stream Mach number where the local Mach number is capable of reaching 
unity. 


Similar results were found for CH-53A conditions. The pressure distribu- 
tions for a high speed level flight condition, case 44, are shown in figure 20 . 
The critical Mach-number effect is apparent in the 30.0 percent chord pressure 
at 85 and 95 percent radial stations but not at 75 percent. Similar to the NH- 
3A case 72, the CH-53A case 44 square wave azimuthal duration is slightly 
greater and begins earlier at r/R = 0.95 than at r/R = O. 85 . This is probably 
due to a higher free-stream Mach number existing toward the blade tip in this 
region. 

Other NH-3A and CH-53A conditions at lower maximum free-stream Mach num- 
bers do not exhibit the step changes in pressure at blade locations that do en- 
counter these changes shown in figures 19 and 20. This also indicates the de- 
pendence on critical Mach number for this effect to occur. The pressures 
reached in these lower Mach-number conditions do not reach those corresponding 
to a critical Mach-number encounter. 

Effects of blade bending. - Reference 6 postulated the possibility of a 
significant effect of flatwise and edgewise blade bending on root torsional mo- 
ments. The use of the product of the flatwise and edgewise bending moments 
measured at a point on the blade, as an indicator gave fairly good phase corre- 
lation of the product and the measured torsional moment retreating blade oscil- 
lations. The product was used as an indicator of load-deflection terms caused 
by the product of blade lift times edgewise blade displacement,; due to bending, 
and by the product of drag load times flatwise bending displacement. No at- 
tempt was made, however, to determine the magnitude of the flatwise/edgewise 
coupling terms. In the present study a semi-empirical method, again utilizing 
the measured flatwise and edgewise moments, was developed to determine the mag- 
nitude, as well as the phase, of the flatwise /edgewise coupling. 

The semi-empirical method developed utilizes a normal mode approach, sim- 
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ilar to that described in reference 8, to model the blade bending. The blade 
was modeled at 15 spanwise segments and was assumed to respond only in its first 
normal mode in the flatwise, edgewise, and torsional degrees of freedom. To 
segregate the flatwise/edgewise coupling influence from other torsional excita- 
tions, the only loads considered to act are the flatwise and edgewise shears 
and moments and the torsional moments at each end of the blade segments . There- 
fore, the torsional moment at the inboard end of any segment can be determined 
from the shears and moments at the, outboard end and the relative flatwise and 
edgewise bending between the ends of the segment. Analytical flatwise and edge- 
wise shears, moments, and displacements can be represented in terms of the blade 
modal response. The measured flatwise and edgewise bending moments on the blade, 
at r/R = 0.65 for the NH-3A and r/R = 0.68 for the CH-53A, were used to describe 
this modal response. From the known boundary conditions at the blade tip, the 
torsion at the inboard end of any segment can be determined from the outboard 
loadings. Therefore, it is possible to derive an equation for the root torsion- 
al moment in terms of blade physical properties, normal mode shapes, and the pro- 
duct of flatwise and edgewise moments measured at one radial station. This semi- 
empirical method was applied to NH-3A and CH-53A cases exhibiting large measured 
retreating blade torsional oscillations. 

Applying this method has shown that, although the phase correlation is 
fairly good, the amplitude of the flatwise/edgewise coupling effect on blade 
torsion is low in the region of the measured oscillations. The comparisons for 
cases exhibiting large measured oscillations are shown in figures 21 and 22. 

The cases shown are high speed NH-3A and CH-53A conditions. For NH-3A case 43 
the maximum coupling oscillation amplitude is only about 20 percent of the cor- 
responding measured torsional oscillation and is less in the regions of the 
other oscillations. The maximum coupling oscillation amplitude for CH-53A case 
1+4 does locally reach approximately 50 percent of the measured value. Other 
coupling oscillations are of lower magnitude however. Additional NH-3A and CH- 
53A conditions exhibit a similar degree of correlation between measured tor- 
sional oscillations and those due to flatwise/edgewise coupling effects. It is 
concluded, therefore, that the coupling of flatwise and edgewise loadings as a 
source of torsional oscillations is small and is not the major cause of measured 
retreating blade torsional oscillations for either of the articulated low twist 
NH-3A (-1+ deg) or CH-53A (-6 deg) blades. In fact, due to the fairly good 
phase correlation, it may be possible for torsion to be a contributor to the 
flatwise and edgewise response of these blades due to coupling. 

Effects of aerodynamic pitching moments. - Discussion in the previous 
sections has ruled out two possible sources as the major cause of retreating 
blade torsional oscillations. The available measured data have shown that no 
one-to-one correspondence appears to exist between possible blade/vortex cross- 
ings and torsional oscillations. Also the coupling of flatwise and edgewise re- 
sponse is not a primary source of the large measured torsional fluctuations. 
Further analysis of the data, however, has shown that the measured oscillations 
appear to be the result of a forced torsional vibration of the blade due to the 
total aerodynamic pitching moment acting on it. This can be demonstrated by ap- 
plying the pitching moment, calculated from the measured blade pressures, to a 
simple single-degree-of- freedom spring/mass torsional model of the blade. As 
an approximation of the aerodynamic pitching moment, the section pitching mo- 



merits calculated from the pressures at 75 , 85 , and 95 percent span stations 
were used to calculate a total moment for the entire blade. A numerical span- 
wise integration was used assuming a linear change in moment from zero at r/R = 
0.30 to the value calculated from the pressures at r/R =0.75. An azimuthal 
time history of the aerodynamic pitching moment was found from this integration. 

A single-degree-of-freedom system with damping was assumed to model the 
blade in torsion. The basic second-order linear differential equation with con- 
stant coeeficients that controls the oscillatory torsional motion at the blade 
root, 0 T , is 


0T + 2c;a>j0.T + to^0 T = j (3) 

where £ is the nondimens ional damping ratio, u>p the torsional natural frequency 
of the blade, T the modal contribution of the aerodynamic pitching moment M x a 
along the blade, and J the torsional modal mass moment of inertia. Only the 
first torsional mode was included in defining uiip , T, and J. A computer program 
was used to numerically solve equation (3) for the root elastic twist 0 rp. Hav- 
ing 0 ip, the model oscillatory root torsional moment M Xj ip was calculated from 

M x t = ^T ®T (^) 

where k^, is the root torsional spring constant. 

The response of the simple torsional model to the applied aerodynamic 
pitching moment was determined for several NH-3A and CH-53A cases. Examples of 
the resulting torsional moment compared to the measured root moment are shown 
in figures 23 and 2b for NH-3A case 71 and CH-53A case UU respectively. The 
curve shown at the top of each figure represents the azimuthal history of the 
aerodynamic pitching moment calculated from the measured blade pressures. Both 
cases indicate the same general trend. A relatively high leading-edge-up moment 
exists in the first quadrant of azimuth, followed by a decrease to approximate- 
ly zero moment on the retreating side and, finally, an increase to the first 
quadrant values. On the advancing side the aerodynamic moment alone matches 
the root structural moment fairly well. The large amplitude torsional oscilla- 
tions experienced by the retreating side structural moment are not exhibited by 
the aerodynamic moment. However, application of the aerodynamic load to the 
model oscillator does result in retreating blade fluctuations of the torsional 
root response. By including 10% of critical damping ( £ = 0.10), it was possible 
to predict fairly well both the magnitude and frequency of the measured oscilla- 
tion. A slight phase lag of the calculated response compared with the measured 
response occurs in both the NH-3A case 71 and the CH- 53 A case UU. This may be 
due to the constant damping factor used in the oscillator model. Actual in- 
flight aerodynamic damping is generally not constant but varies with blade 
azimuth. 

The mechanism causing the retreating blade torsional moments is believed 
to be initiated by the advancing blade aerodynamic pitching moment. The sudden 
decrease in large leading-edge-up advancing blade moment appears to initiate 
oscillations at the blade torsional natural frequency for the remainder of the 
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rotor cycle. These occur due to simple blade spring and inertia effects. The 
NH-3A first torsional natural frequency is 6.6/rev while the CH-53A's is 6.0/rev. 
Variations of the aerodynamic moment later in the cycle either tend to increase 
the amplitude of the oscillations or to assist in damping them out. Although 
the calculation of the aerodynamic pitching moment was approximate and a simpli- 
fied single-degree-of-freedom blade model was used, the prediction of oscilla- 
tions due to the application of only the aerodynamic moment is encouraging. 


Correlation with Analytically Derived Results 

The purpose of the correlation study was both to evaluate the present an- 
alytical capability for predicting helicopter blade torsional response and to 
understand better the nature of the torsional excitation for maneuvering and 
high-speed flight conditions. The effects of unsteady aerodynamics and non-uni- 
form rotor inflow upon the correlation of blade section aerodynamic coefficients 
and blade structural torsional moment were evaluated. The analytical solutions 
were examined to determine the sources of the calculated torsional excitation. 

Description of programs and procedures .- Analytical description of rotor 
blade response requires proper modeling of both dynamic and aerodynamic effects. 
The Normal Modes Aeroelastic Blade Analysis was used to calculate blade elastic 
response. As described in reference 8 this analysis represents blade flatwise, 
edgewise, and torsional elastic deformation by a summation of normal mode re- 
sponses. For a given rotor inflow and a set of control positions this program 
carries out a timewise integration of the blade equations of motion to solve for 
the modal response. Blade elastic motions and aerodynamic forces and moments 
are fully coupled. At each radial and azimuthal position, blade section loading 
and pitching moment coefficients are obtained from appropriate steady or unsteady 
data tables. To model unsteady aerodynamic effects the procedures and data 
given in reference 6 were used. Briefly this method is based on airfoil data 
obtained from a 0012 airfoil section undergoing sinusoidal angle-of-attack 
changes. Section lift and pitching moment coefficients were measured and tabu- 
lated as functions of instantaneous blade angle of attack and its first two time 
derivatives. Unsteady aerodynamic lift and moment coefficients were obtained in 
the computer program by entering the unsteady data tabulations with calculated 
angle-of-attack parameters. Because only unsteady 0012 data are available, they 
were used to approximate the unsteady 0011 section characteristics of the CH-53A. 
Steady drag data were used throughout the study. Variable rotor inflow was de- 
termined using the Prescribed Wake Analysis which is described in reference 9» 

For a given flight condition this program calculates the induced velocity distri- 
bution over the rotor disc due to the rotor wake. For this study the wake was 
approximated by the classical skewed helicoid defined when the trailing vortex 
filaments generated by the blade are carried downstream. When the positions 
and strengths of the vortex filaments are determined the velocity induced at 
each blade position is calculated using the Biot-Savart law. In this way close 
blade-vortex passages result in large and rapid fluctuations in induced velocity. 
The program is strictly a steady-flight analysis. Its use for maneuvers implies 
a quasi-steady treatment of the time varying aircraft cqndition. 

In the aeroelastic analysis, blade flexibility was modeled using rigid 
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body flapping and lagging degrees of freedom and six elastic modes: three flat- 

wise, two edgewise, and one torsional. For the level-flight conditions, use of 
measured cyclic and collective pitch and rotor shaft angle gave analytical rotor 
lift and propulsive force very near test values. In maneuvering flight, where 
rotor shaft angle was not accurately known, an iteration was performed in which 
shaft angle was varied to achieve approximately the correct rotor lift. Since 
both unsteady aerodynamics and variable rotor inflow were thought to be necessary 
for modeling high-frequency rotor torsional response, the separate effects of 
unsteady aerodynamics and variable inflow were examined. Accordingly the aero- 
elastic program was run with steady aerodynamics and variable inflow, unsteady 
aerodynamics and variable inflow, and with unsteady aerodynamics and constant in- 
flow corresponding to the average momentum induced velocity. 

In order to establish representative trends for the effects of aerodynamic 
and inflow modeling and to assess the degree of correlation, six cases were 
studied in detail. Of the conditions shown in table I, NH-3A cases 43 and 72 
and CH-53A cases 4lLF, 4lTR, 4lTN, and 44 were analyzed and the resulting section 
loading and pitching moment coefficients, total aerodynamic blade pitching mo- 
ments, and root torsional moments were plotted and compared with test data. It 
was felt that the results of NH-3A case 72 and CH-53A cases 4lLF, 4lTN, and 44 
were representative and the results presented in the following section will be 
taken entirely from these cases. 

Discussion of results .- The results of the correlation study indicate 
that rotor blade aerodynamic /dynamic behavior can be reasonably well predicted 
using the present analytical methods. Certain problem areas are pointed out 
however. Insight into the mechanism of torsional motion is also provided. 

For NH-3A case 72, which is a 60 deg angle of bank turn at m = 0.35 and 
C t/o = 0.085, a full set of correlation results is given in figures 25 through 
28. The abbreviations V.I. and C.I. used in figures 25 through 42 refer respec- 
tively to the use of a variable inflow and a constant inflow analysis. In figure 
25, which shows the root torsional moment correlation, three results are evident. 

First the test data show a steady positive displacement from the analytical re- 

sults of approximately 400 N-m (3540 in. -lb). Second the use of variable inflow 
in place of the uniform momentum velocity greatly improves the correlation of 
the high-frequency response. Finally the use of unsteady aerodynamics led to 
slightly higher retreating blade torsional oscillations although test values 
were not achieved. Examination of figure 26 which shows the correlation of total 
analytical aerodynamic pitching moment and pitching moment calculated from test 
pressures indicates that the test aerodynamic moment is more positive than the 

analytical predictions. Also the higher frequency makeup of the time history is 

best matched when variable rotor inflow is used. The buildup of torsional mo- 
ment for the variable inflow cases appears to be the result of local blade stall 
which begins inboard at r/R = 0.50 at ip = 180 deg and propagates out along the • 
blade through the third quadrant of rotor azimuth. Figure 27 shows the compari- 
son of the outboard c m coefficients reduced from test data and those calculated 
in the Normal Modes Analysis. As described previously the noseup pitchingjno- 
ment coefficients which are evident in the test data are believed to be the result 
of asymmetry in the airfoil contour introduced by pressure transducer fairing. 

The point of rapid pitching moment coefficient buildup is fairly well matched by 
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the analytical results, although the c m variation in stall is not well described. 
Figure 28 shows the time histories of blade section loading coefficients at r/R 
= 0.75, 0.85, and 0.95. The best correlation with test data is achieved using 
unsteady aerodynamics and variable rotor inflow. 

CH-53A case 4l was studied during the level-flight, transition, and steady- 
state turn phases. These data were taken at y = 0.24 and Cr p/a ranging from 0.063 
for level flight to 0.099 for the turn. Results of the level flight and steady- 
state turn calculations are shown in figures 29 through 32 and 33 through 36 re- 
spectively. The two root torsion moment plots, figures 29 and 33, show that the 
rise in oscillatory control moment is predicted by the analysis. As observed 
for the NH-3A case the measured moments are generally more noseup. In the steady- 
state turn the rotor angle of attack relative to the free stream was not known 
and a shaft angle of attack of 0 deg was used to make the analytical lift agree 
with the measured value. One result of the small shaft angle was that the rotor 
wake passed through the cone prescribed by the blades, causing large induced 
velocity spikes at inboard blade sections and precipitating local stall an d 
large nosedown aerodynamic pitching moments. Figures 30 and 34 show the dif- 
ferences between the aerodynamic pitching moment calculated analytically and 
that determined by the approximate integration of section pitching moment data. 

The integrated test aerodynamic moment for this case using only outboard pitch- 
ing moment data does not agree well with the analytical calculation. Samples of 
the section pitching moment coefficient correlation for the level-flight and 
turn conditions are shown in figures 31 and 35. As previously shown for the NH- 
3A the test pitching moment data show a steady noseup value which is not pre- 
dicted by the analysis with standard airfoil data. Typical section loading co- 
efficient results shown in figures 32 and 36 depict the increase in the vibratory 
component of blade lift which accompanies higher blade loading and deeper pene- 
tration into stall. For case 4lTN , the variable inflow program has apparently 
overpredicted the level of velocity fluctuation felt by the advancing blade. In- 
accurate modeling of either the strength or position of the trailed vortices is 
the apparent source of the discrepancy. 

The final correlation results to be presented are for CH-53A case 44 which 
was a level-flight condition flown at y = 0.38 and C T /a = 0.057. Figure 37 shows 
that none of the analytical results predict the high level of test oscillatory 
root torsional moment. For this case it was shown in figure 24 that the blade 
root torsional signature could be fairly well described as the forced response to 
the applied aerodynamic pitching moment. The analytical pitching moment coeffi- 
cients shown in figures 38 and 39 &o not match the test values. The solutions 
obtained using unsteady aerodynamics and variable inflow generally give the best 
agreement but do not match the test advancing blade positive pitching moment co- 
efficients. Fairly good agreement between test and analytical normal force coef- 
ficients is shown in figures 40 and 4l. 

As a further investigation of the torsional response for CH-53A case 44, 
an attempt was made to force the analytical aerodynamic moment to agree with that 
calculated from test pressures. In this way, calculated blade dynamic response 
could be compared directly with measured response for the same input moment. Ac- 
cordingly the standard 0011 c m data were modified by shifting the below-stall 
pitching moment coefficients to 0.02. This value was chosen as an average of 
the c m values in the test data. The changes in aerodynamic pitching moment and 
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root torsional moment caused "by this modification are shown in figure 1*2. Al- 
though the desired nose-up change in advancing blade aerodynamic moment was only 
partly accomplished, an improvement in both the amplitude and azimuthal variation 
of the resulting torsional moment is achieved. This supports what was found 
above with the simple single-degree-of-freedom torsional model. Torsional mo- 
ments applied to the advancing blade serve as initial conditions for retreating 
blade response. Without proper modeling of advancing blade aerodynamics the re- 
treating blade torsional oscillation cannot be obtained. 

A final attempt to understand the nature of the blade torsional response 
was made by examining the terms which contribute to analytical torsional response. 
In the Normal Modes Analysis coupled equations for blade flatwise, edgewise, and 
torsional modal acceleration are solved simultaneously at each rotor azimuth. 

For the CH-53A case kk, run with the modified pitching moment data, the time his- 
tories of the major terms contributing to non-dimensional torsional modal re- 
sponse are shown in figure U3. The three primary effects which determine the re- 
sponse of the first elastic torsional mode are (l) the aerodynamic contribution 
which is made up of pure pitching moment and moments created around the blade 
elastic axis by lift and drag acting at the blade quarter chord, (2) the blade 
torsional inertia, and (3) dynamic moments about the blade elastic axis caused 
by flapping and centrifugal acceleration acting at the blade section centers of 
gravity. The high-frequency torsional modal response in the fourth quadrant is 
primarily the result of the applied aerodynamic moment and the blade torsional 
inertia. In agreement with the results shown above for the simple torsional os- 
cillator model, this shows that the blade torsional moment, which is proportional 
to the modal response, Is basically a forced response to the applied aerodynamic 
moment . 


ANALYSIS OF ROTOR VIBRATORY FORCE DATA 


A secondary area of study included in this program involved an analysis of 
the in-flight vibratory loads encountered by the main rotor on both the NH-3A 
and the CH-53A, These loads were calculated from measurements recorded during 
earlier flight-test programs over a wide range of operating conditions. The 
study reported herein included a determination of the primary dynamic terms re- 
sulting in the harmonics of blade hinge forces and the b/rev hub loads. Also 
the effects of the number of main rotor blades on harmonics of the hinge forces 
were sought. 


Major Blade Hinge Force Dynamic Terms 

The purpose of this study was to attempt to define, in a relatively few 
number of terms, a method to calculate the harmonics of main rotor articulated 
blade hinge forces from blade geometry and test measurements. Knowledge of the 
harmonic content of the hinge forces enables the determination of the vibratory 
loads fed by the rotor into the fuselage. Reference h describes the successful 
use of the Root Response Method to determine the. hinge forces on the NH-3A over 
a wide range of flight conditions. This method was also used successfully with 
the CH-53A flight data. To determine the loads, measurements near the root end 
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of one blade are recorded during flight. From these data a computer program 
calculates the azimuthal history and corresponding Fourier series coefficients 
of the three rotating system hinge forces: radial, F£; edgewise, F^; and ver- 

tical, F*. Since measurements are made in the local blade axis system, trans- 
formation of the equations of motion between the local blade and rotor head 
axis systems is required. The axis systems involved are shown in figure 44 . 

The resulting hinge force equations in the rotor head system include between 80 
and 150 dynamic terms. Definition of the terms are included in reference 4 . It 
is desirable to simplify these expressions to include only the major terms. 

During the present study, analysis of the available NH-3A and CH-53A 
flight data has resulted in a simplified expression for each of the three hinge 
forces. A maximum of only six terms was found necessary to describe any of the 
forces on either aircraft. The number of terms required to give an approximation 
of any harmonic from the first through the seventh, generally within 15 percent 
in amplitude and 10 deg in phase, varies between one and six. Table II indicates 
the terms required to achieve this accuracy for NH-3A and CH-53A rotor systems. 
The first seven harmonics are included since they comprise all harmonics required 
to determine the b/rev hub loads in the fixed system of either aircraft. Gen- 
erally the same terms, with the exception of the X3, Y5, Y7, and Z 6 hub motion 
terms, are important for the same harmonic on both aircraft. The hub motion 
terms are important for different harmonics due to the difference in the number 
of main rotor blades, i.e., five on the NH-3A and six on the CH-53A. As table 
II indicates, by measuring the flatwise and edgewise bending at each end of a 
short root end instrumented segment, the blade tension near the root, the root 
blade angular motions (hunting, flapping, and feathering), and the root head 
translation motions, the harmonics of the hinge forces can be determined. 

Tables III and IV indicate the application of these approximations to 
typical NH-3A and CH-53A conditions. Compared are the summations of the major 
terms, indicated in table II, and the exact calculation comprised of all terms. 
Harmonic error of the approximation is generally within 15 percent of the exact 
amplitude and 10 deg of the correct phase. Also apparent in each direction is 
the general decrease in harmonic amplitude with increased harmonic order. 

An example of how only a few of the many terms contributing to the calcu- 
lation of the hinge force can approximate the exact calculation is shown in 
figure 45. Shown are the major terms contributing to the fifth harmonic of the 
edgewise hinge force for a high-speed CH-53A condition. The five terms required 
to approximate the harmonic, indicated in table II, are plotted in vector form 
along with their vector sum, F^- 5 j^, and the exact calculation, F^ 5 . In this 
instance the approximation is within 7 percent of amplitude and 5 &eg of the 
exact phase. Each of the five terms is of significant amplitude. The source of 
the Y1 term is the product of measured blade tension and root lag motion. The 
edgewise shear force causing edgewise blade bending across the root instrumented 
segment results in the largest magnitude term, Y2. The Y 6 and Y 8 terms exist 
due to lag motion inertia of the blade cuff mass inboard of the instrumented 
segment. The source of the Y7 term is the vibratory motion of the rotor head. 
Other harmonics indicate similar correlation between the exact and approximate 
solutions . 
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It has been shown that, for an articulated rotor system, the vibratory 
loads acting on the rotor head in the rotating system can be determined fairly 
accurately using only a few dynamic terms. Azimuthal histories of measurements 
near the root of one blade can be recorded and the terms indicated in table II 
calculated in time history form. Standard harmonic analysis of these time his- 
tories then provides the required Fourier coefficients of each term. 


Major Rotor Huh Force Dynamic Terms 


Knowledge of the blade hinge forces encountered in flight enables one to 
calculate the fixed system forces which the rotor system feeds into the fuselage 
and which consequently excite the airframe . The fixed system vibratory loads 
occur, when each blade is assumed to have the same azimuthal hinge forces as 
does the instrumented blade, at a frequency that is an integer multiple of the 
number of blades times the rotor speed. A 5-bladed rotor system for example, 
can theoretically excite the fuselage with loads at 5/rev, 10 /rev, 15/rev, etc. 
For articulated rotors, generally only the b/rev loads have significant ampli- 
tudes compared to the 2b/rev, 3b/rev, etc. loads. Therefore, knowing the major 
dynamic terms contributing to the rotating system hinge force harmonics, an 
analysis was undertaken to determine which of these were the primary contribu- 
tors to the fixed system b/rev hub loads. Specifically the major dynamic blade 
terms contributing to the b/rev hub forces were defined for both the NH-3A and 
the CH-53A rotor systems. 


The b/rev hub forces, Fx,bi Fy^bj and Kz,b can represented as Fourier 
cosine and sine components of the b-th harmonic in terms of the hinge force 
harmonic coefficients as follows 


F = — (f* + F* + F* - F* \ 

X,bc 2 x,b-l,c x,b+l,c y,b-l,s y,b+l,s) 

F = — (f* + F* _ F* + F* % 

X, bs 2 x,b-l,s x,b+l,s y,b-l,c y,b+l,c) 

F = — (-F* +F* + F* + F* \ 

Y, bc 2 x,b-l,s x,b+l,s y,b-l,c y,b+l,c) 

F = — (F* - F* + F* + F* \ 

Y,bs 2 v x,b-l,c x,b+l,c y,b-l,s y,b+l,s) 


(5) 

( 6 ) 
(T) 

( 8 ) 


F = bF* . 

Z,bc z,b,c 

F = bF* 

Z,bs z,b,s 


(9) 

( 10 ) 


where the hinge force cosine and sine harmonic coefficients are defined from 
the corresponding amplitude and phase. For example 


F * V 

z,b,c 

= f# 
z,b 

sin<}> 

F* , 

= F* 

cos<}> 

z,b,s 

z,b 


( 11 ) 

( 12 ) 
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The same NH-3A and CH-53A flight conditions that were used to determine the 
major hinge force dynamic terms were also analyzed to determine which of these 
terms are also the primary terms causing the fixed system b/rev hub forces. As 
equations 5 through 10 indicate, only the (b±l)-th in-plane hinge forces (F* and 
F^) and the b-th vertical hinge force (F*) excite the fixed system at the b/rev 
frequency. It was found that generally all of the terms required to achieve a 
good approximation of the hinge force harmonics are also needed to give a 
reasonable estimation of the b/rev hub forces. 

Figures 46 through 48 and 49 through 51 indicate, respectively for NH-3A 
and CH-53A conditions at y = 0.36, the major dynamic terms contributing to the 
hub forces. Also included is the degree of agreement with the exact calculation. 
The upper chart in each figure indicates the major terms contributing to the 
b/rev cosine component of the fixed system force. The lower chart shows those 
for the sine component. The NH-3A 5/rev in-plane hub forces result from 4 and 
6/rev rotating system in-plane hinge forces while the 5/rev vertical hub force 
comes from 5/rev vertical hinge forces. Similarly, 5 and 7/rev in-plane and 
6/rev vertical hinge forces transform into the 6/rev CH-53A fixed system loads. 
Comparison of the major NH-3A and CH-53A terms shows in many cases that there 
are usually several large magnitude terms causing the fixed system vibratory 
loads. Also the same terms do not always exhibit the largest magnitudes and, . 
therefore, all terms shown in figures 46 through 51 must be included to form a 
reasonably good approximation of the exact calculation. For the NH-3A example 
shown the approximations indicated give excellent agreement with the exact com- 
putation. The CH-53A case indicates an acceptable correlation, but not as good 
as the NH-3A example for some forces. This represents the minimum degree of 
correlation between the exact and approximate computations for the CH-53A condi- 
tions analyzed. The vertical loads in figure 51 indicate the poorest degree of 
correlation. However, comparison with the corresponding NH-3A loads in figure 
48 and other CH-53A conditions indicates that the magnitude of these loads is 
low. 


Effects of the Number of Main Rotor Blades on 
Harmonics of Blade Hinge Forces 

The availability of blade hinge force data for both the 5-bladed NH-3A 
and the 6-bladed CH-53A over a wide range of operating conditions made possible 
a study of the effects of the number of blades on these loads. Available har- 
monics of the hinge forces extended from level-flight conditions at y = 0.09 to 
y = 0.36 for both aircraft. The flight conditions used for the NH-3A correspond 
to a nominal gross weight of 8560 kg (18900 lb) and a neutral aircraft center 
of gravity. The NH-3A was in a compound configuration (auxiliary jets and wing 
installed) but the conditions to be discussed were flown with the jets at idle. 
A nominal gross weight of 15860 kg (35000 lb) and neutral aircraft center of 
gravity comprised the CH-53A configuration. Comparisons of the first seven 
harmonics of the NH-3A and CH-53A hinge forces are shown in figures 52, 53, and 
54 at advance ratios of 0.1, 0.2, and 0.3. Each force harmonic has been non- 
dimensionalized by the steady (mean) value of vertical hinge force. 

Figures 52, 53, and 54 show that generally higher amplitude harmonics 
exist for the CH-53A with the primary exception being the first harmonic verti- 
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cal force. This is due to differences in rotor trim for the conditions con- 
sidered. Similar trends occur at the three advance ratios shown. The two in- 
plane hinge forces, F£ jn and Fy n ,. generally have relatively low amplitude 
second harmonics, higher middle order harmonics, and low amplitude high order 
harmonics (sixth and seventh). The harmonic amplitude of the vertical hinge 
force generally decreases as harmonic order increases . The large amplitude 
fourth harmonic edgewise hinge force on the CH-53A, shown in figure 53, is due 
to an edgewise blade resonance near the corresponding k/vev frequency. 

It is believed that figures 52, 53, and 5k can be used to approximate the 
magnitude of hinge force harmonics generated by an articulated rotor system if 
they are known for another hinged rotor system having one more or one less blade. 
This should be applicable to blades of similar construction and geometry as the 
NH-3A or CH-53A blades. Applicability of these curves to blades having consid- 
erably different construction, geometry, or number of blades per rotor requires 
further investigation. 


CONCLUSIONS 


1. As anticipated, larger high-frequency control loads were found to occur 
on the retreating blade as rotor load factor or airspeed increased. 

2. From the test data analyzed, only a limited blade/tip vortex interaction 
effect appears in the measured blade pressures and resulting section load- 
ing and pitching moment coefficients. The CH-53A blade/vortex interactions 
appeared to have a stronger effect than those on the NH-3A. The most no- 
ticeable effects were experienced by each aircraft in maneuvering, rather 
than level flight, conditions. Even here, however, these interactions 

did not appear to dominate the torsional moment oscillatory amplitudes. 

3. Blade/tip vortex interactions do not exhibit a direct azimuthal correla- 
tion with retreating blade torsional oscillations. 

4. Evidence of blade upper surface shock, due to a critical Mach-number en- 
counter, was found in the measured pressure data, corroborating similar 
observations by NASA on the H-3*+ blade. 

5. Coupling of flatwise and edgewise blade response is not a major cause of 
the measured retreating blade torsional oscillations for the low twist 
blades tested and analyzed. 

6. The high-frequency blade control loads primarily result from the forced 
dynamic response of the blades produced by the time varying aerodynamic 
pitching moments. Although fluctuations in the pitching moment caused by 
stall flutter effects may serve to amplify or attenuate retreating blade 
response, the relatively rapid changes in pitching moments on the advanc- 
ing blades appear to excite blade torsional oscillations that persist to 
the retreating side of the disc. (Stall flutter is probably a major con- 
tributor to high-frequency control load oscillations when operating well 
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above the knee of the control load curve, but it was found that for the 
less stalled conditions analyzed in this study, the principal cause of 
high-frequency control loads appears to be a forced dynamic response of 
the blade due to the aerodynamic moment.) 

7. The instrumented blades exhibited significant positive pitching moments 
in flight test despite the fact that the airfoil section was nominally 
symmetrical. These moments may have been the result of instrumentation 
fairings used or possibly may have arisen from a more fundamental source. 

8. Fair correlation of measured blade airloadings with those predicted using 
a normal mode, rigid hub, blade aeroelastic analysis was achieved. The 
inclusion of unsteady aerodynamics and variable inflow effects results in 
the best agreement. To improve the correlation it appears mandatory to 
include in the analysis the positive airfoil pitching moment coefficients 
observed in the test data.. 

9. The equations presented in reference 4 for computing hub vibratory forces 
from measurements of blade root dynamic response can be considerably sim- 
plified. Many of the dynamic correction terms due to motion of the hub, 
blade cuff, and instrumented blade segment can be neglected. 

10. For the two rotors analyzed, the harmonics of root shear, nondimens ional- 
ized by steady axial (vertical) blade shear, differed significantly, in- 
dicating the need for (l) further work to define all controlling parameters 
and (2) caution in applying these results to other rotor designs. 

RECOMMENDATIONS 


1. ' Additional work should be conducted to investigate the source of the 

positive pitching moments observed on the flight test blades. At the 
very least, two dimensional tests of a simulated blade section should be 
conducted to determine if the instrumentation fairings were a significant 
contributor. 

2. Improved procedures for modeling the airloads produced on a blade passing 
close to a series of vortices should be developed based on a lifting sur- 
face representation of the blade. 

3. A study should be conducted to correlate vibratory root shears predicted 
by analysis, with those measured. 


Sikorsky Aircraft Division, 

United Aircraft Corporation, 

Stratford, Connecticut, November 21, 1972 
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TABLE I. - TEST CONDITIONS FOR HIGH-FREQUENCY BLADE RESPONSE STUDY 
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DYNAMIC TERMS REQUIRED FOR APPROXIMATION OF NH-3A AND 
CH-53A HARMONICS OF HINGE FORCES 
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Terms Required for Approximation of Total Harmonic of Hinge Force 
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TABLE II. - DYNAMIC TEEMS REQUIRED FOR APPROXIMATION OF NH-3A AND 
CH-53A HARMONICS OF HINGE FORCES - Concluded 
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TABLE III . - SAMPLE COMPARISON OF APPROXIMATE AND EXACT CALCULATION 
OF NH-3A HARMONICS OF BLADE HINGE FORCES 


Compound mode ,V =72.0 m/sec (l40;0 kn), y 


Harmonic of 
Hinge Force 


0.36, C T /o = 0.0lt3 




Summation of Major 
Dynamic Terms 
(Approximate Solution) 

Summation of all Dynamic 
Terms 

(Exact Solution) 

Amplitude 

Phase 

Amplitude 

Phase 

N 

lb 

deg 

N 

lb 

deg 

2 1*80 

559 


2520 

568 

228 

1*31 

97 

■Si ■ 

534 ! 

120 

275 

516 

116 

l4o 

516 

116 

135 

1000 

226 

349 

1050 

237 

354 

271 

6l 

1 66 

298 

67 

163 

725 

163 

308 

725 

163 

307 

85 

19 

115 

95 | 

21 

114 

1940 

436 

58 

1930 

433 

62 

71*7 

168 

13 

770 

173 

13 

587 

132 

348 

578 

130 

350 

551 

124 

271 

525 

118 

271 

56 

13 

1 

50 

11 

11 

ll+2 

32 

21 

162 

37 

27 

31* 

8 

355 

36 

8 

359 

231*00 

5270 

103 

24400 

5470 

102 

11*60 

328 

253 

1420 

320 

249 

1990 

448 

211 

2080 

467 

214 

83 

19 

7 

102 

23 

359 

322 

72 

209 

300 

67 

210 

316 

71 

327 

313 

70 

326 

34 

8 

154 

34 

8 

150 


Harmonic of Hinge Force 
Absolute Error Due to 
Approximation 
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TABLE IV. - SAMPLE COMPARISON OP APPROXIMATE AND EXACT CALCULATION 
OF CH-53A HARMONICS OF BLADE HINGE FORCES 


V = 78.4 m/ sec (152.4 kn), p = 0.36, C^/a = 0.055 


Harmonic of 
Hinge Force 

Summation of Major 
Dynamic Terms 
(Approximate Solution 

Summation of All Dynamic 
Terms 

(Exact Solution) 

Harmonic of Hinge Force 
Absolute Error Due to 
Approximat i on 


Amplitude 

Phase 

Amplitude 

Phase 

Amplitude 

Phase 


N 

lb 

deg 

N 

lb 

deg 

% 

deg 

F* 

X,1 

6510 

1470 

182 

7160 

1610 

182 

9 

0 

F* 

x ,2 

582 

131 

129 

595 

134 

128 

2 

1 

F* , 
x,3 

1510 

340 

356 

1740 

391 

353 

13 

3 

F x,4 

1900 

427 

340 

2090 

470 

332 

9 

8 

F x,5 

831 

187 

248 

685 

154 

235 

21 

13 

F* C 
x ,6 

818 

184 

354 

916 

206 

1 

11 

7 

F* 

x,T 

268 

60 

324 

266 

60 

318 

1 

6 

F* 

y»i 

4400 

990 

67 

4850 

1090 

70 

9 

3 

F* 

y .2 

480 

108 

153 

511 

115 

154 

6 

1 

F* 

y>3 

2010 

452 

282 

2350 

528 

278 

14 

4 

y>4 

3120 

700 

256 

3460 

779 

256 

10 

0 

F* „ 
y»5 

408 

92 

59 

44o 

99 

54 

7 

5 

F* A 
y .6 

595 

134 

32 

734 

165 

25 

19 

7 

F* 

y»T 

249 

56 

79 

295 

66 

74 

16 

5 

F* 
z ,1 

16600 

3720 

92 

17200 

3870 

88 

4 

4 

F* 

z ,2 

6750 

1520 

253 

6540 

1470 

251 

3 

2 

F* _ 
z,3 

7250 

1630 

155 

7470 

1680 

159 

3 

4 

F* , 
Z ,4 

911 

205 

353 

782 

176 

344 

16 

9 

F* _ 
Z,5 

1160 

260 

250 

1300 

293 

253 

11 

3 

F * (L 
z ,6 

130 

29 

173 

101 

23 

206 

29 

33 

F* 

z,7 

149 

34 

26 

176 

40 

9 

15 

17 
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Figure 1.- NH-3A Test Aircraft. 
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Figure h.- Blade Pressure Transducer Locations 
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Figure 8— Maneuver Effects on CH-53A Main Rotor Blade Section Pitching 
■ Moment Coefficient at r/R = 0.85. 
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Figure 9-,~ Maneuver Effects on CH-53A Root Torsional Moment 
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Figure 10. _ Plan View of Non-distorted Wake Helical Pattern for NH-3A Case 72, Right Turn 
60 deg Angle of Bank, y = 0.35. 



Figure 11.- Azimuthal Histories of Blade Differential Pressure at r/R 
for NH-3A Case 72. 







igure 12.- Azimuthal Histories of Blade Differential Pressure at r/R =0.85 
for NH-3A Case 72. 
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Figure 15 •- Azimuthal Histories of Blade Differential Pressure at r/R = 0.75 
for CH-53A Case 4lTN . 






o 



>.01 X 2 uj/n ‘ d V ‘ 3dnSS3dd lVHN3a3dJIQ 


■ ■ i i 1 

ID CM 00 ^ O * 


LTN 

co 

o 

II 


ISd‘dV‘3dnSS3dd lVHN3d3ddia 


Figure 16 .- Azimuthal Histories of Blade Differential Pressure at r/R 
for CH-53A Case klTN . 
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Figure IT— Azimuthal Histories of Blade Differential Pressure at r/R = 0.95 
for CH-53A Case UlTH. 





ROOT TORSIONAL MOMENT , M x , IN.- LB 


24 OOO r 


HH 


ISWfll 


POSSIBLE VORTEX 
Ci CROSSINGS (T) 
WITH r/R =75 


ptm 


uni 


nil 


aBHW 


AZIMUTH , , DEG 


Figure l8.- Azimuthal Histories of Section Loading Coefficient, Section 
Pitching Moment Coefficient, and Root Torsional Moment for 
CH-53A Case 4lTN. 
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Figure 22.- Contribution of Blade Flatwise and Edgewise Coupling to Root 
Torsion for CH-53A Case U4. 
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Figure 23.- Comparison Between Measured Blade Torsional Response and Torsional Oscillator Model 
Response to Aerodynamic Pitching Moment for NH-3A Case 71. 
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Comparison Between Measured Blade Torsional Response and Torsional 
Oscillator Model Response to Aerodynamic Pitching Moment for CH-53A 
Case 4U. 
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Figure 25-- Correlation of Root Torsional Moment for NH-3A Case 72 
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Figure- 26.- Correlation of Aerodynamic Pitching Moment for NH-3A Case 72 
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(b) r/R =0.85 


Figure 28.- Continued. 
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29.- Correlation of Root Torsional Moment for CH-53A Case UlLF, 
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Figure 30.- Correlation of Aerodynamic Pitching Moment for CH-53A Case HlLF. 
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Figure 31.- Correlation of Section Pitching Moment Coefficient at r/R 
for CH-53A Case 1+1LF. 
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Figure 33.- Correlation of Root Torsional Moment for CH-53A Case 4lTN 
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Figure 3^.- Correlation of Aerodynamic Pitching Moment for CH-53A Case 






Figure 35 . - Correlation of Section Pitching Moment Coefficient at r/R 
for CH-53A Case 1*1 TN. 
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Figure 36.- Correlation of Section Loading Coefficient at r/R = O .85 for 
CH-53A Case 4lTU. 
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Figure 37*- Correlation of Root Torsional Moment for CH-53A Case U4 
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Figure 38 .- Correlation of Section Pitching Moment Coefficient at r/R = 0.75 
for CH-53A Case 1+4. 
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Figure 39*- Correlation of Section Pitching Moment Coefficient at r/R = 0.95 
for CH-53A Case ^4. 
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Figure 40. — Correlation of Section Loading Coefficient at r/R = 0. 75 for 
CH-53A Case 44. 
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Figure 4l.- Correlation of Section Loading Coefficient at r/R = 0.95 for 
CH-53A Case 1*U. 
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Figure ^2.- Effects of Modifications to Analytical c Aerodynamic Data on 

Aerodynamic Pitching Moment and Root Torsional Moment for CH^53A 
Case UH. 
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Figure 43.,- Major Analytical Contributions to Nondimensional Blade. Torsional 

Response for CH-53A Case 44 With Unsteady Aerodynamics, Variable 

Inflow, and Modified c Data. 
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COMPLETE DEFINITIONS OF THE 
SYMBOLIC TERMS <Y‘s ) ARE 
GIVEN IN TABLE II . 



Figure ^5-- Comparison of Approximate and Exact Calculations of the Fifth 
Harmonic of Edgewise Hinge Force, CH-53A, V = J8.h m/ sec 

(152. U kn), y = 0.36, C T = 0.055- 
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Figure 51 •- Major b/Rev Vertical Hub Force Dynamic Terms, CH-53A 
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Figure 53 .- Comparison of NH-3A and CH-53A Harmonics of Edgewise Hinge Force. 
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